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Typical density profile for warm dark matter haloes 
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ABSTRACT 

Using the model recently developed by Salvador-Sol et al. (2012), we derive the typi- 
cal spherically averaged halo density profile from the power-spectrum of density per- 
turbations in the concordant A warm dark matter (WDM) cosmology with 2 keV 
non-thermal sterile neutrinos. This allows us to analyse separately the effects on the 
density profile at small radii of the spectrum cutoff caused by free-streaming and the 
bound in the fine-grained phase-space density, both due to the non-negligible particle 
velocities at decoupling. 
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1 INTRODUCTION 

Matter in the Universe is predominantly dark and clustered 
in haloes that grow through mergers and accretion. The 
concordant A, dissipationless collisionless, cold dark matter 
(CDM) model recovers the observed large-scale properties 
of the Universe: it correctly predicts the microwave back- 
ground radiation anisotropies (Komatsu et al. 2011) and 
galaxy clustering (Cole et al. 2005). However, some problems 
arise in the small-scale regime: it predicts excessive substruc- 
ture with a deficient distribution of maximum circular veloc- 
ities (Klypin et al. 1999; Moore et al. 1999; Boylan-Kolchin 
et al. 2011) and a sharp central cusp in the halo density pro- 
file, apparently in conflict with the profiles observed in dwarf 
galaxies (Goerdt et al. 2006). While the disagreement in the 
satellite abundance and characteristics might be explained 
through the inhibition of star formation owing to several 
baryonic feedback processes, the cusp problem seems to be 
insurmountable. 

Several modifications of the ACDM model have been 
proposed that, keeping its right predictions at large 
scales, may improve those compromised at small scales. 
This includes the cosmological models dominated by self- 
interacting dark matter (Spergel & Steinhart 2000; Kapling- 
hat et al. 2000; Bento et al. 2000) and dissipationless col- 
lisionless warm dark matter (WDM). The best candidate 
particles in the latter category of dark matter are the grav- 
itino (Ellis et al. 1984; Hogan & Dalcanton 2000; Kaplinghat 
et al. 2005; Gorbunov et al. 2008) and non-thermal sterile 
neutrino (Dodelson & Widrow 1994; Hogan & Dalcanton 
2000; Shaposhnikov & Tkachev 2006). Specifically, the case 
of light sterile neutrinos is being attracting growing interest 
as this kind of particles is naturally foreseen in a minimal 
extension of the Standard Model. 
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The non- negligible WDM particle velocities at decou- 
pling introduce a cutoff in the power-spectrum due to free- 
streaming and, at the same time, a bound in the fine-grained 
phase-space density. The former effect should inhibit the 
formation of haloes below the corresponding free-streaming 
mass scale, Mf a , which would presumably produce a cen- 
tral core in the subsequent generations of more massive 
haloes. The latter should translate into an upper bound in 
the coarse-grained phase-space density of haloes resulting 
from virialisation, which could also lead to the formation 
of a core. Both aspects depend, however, critically on the 
poorly known way haloes fix their density profiles. 

In fact, A-body simulations do not confirm such expec- 
tations: WDM haloes show similar density profiles as CDM 
haloes (e.g. Colombi et al. 1996; Colin et al. 2000; Knebe et 
al. 2002; Wang & White 2009; Schneider et al. 2011). Only 
a small hint towards increased scaled radii has been found 
(Bode at al. 2001), although the opposite trend, namely 
a slight inflection towards sharper density profile at small 
radii, has also been reported (Colin et al. 2008). The sit- 
uation is further complicated by the fact that simulations 
of WDM cosmologies find a substantial amount of haloes 
with masses considerably smaller than M{ B . These low-mass 
haloes may be spurious due to the effects of the periodical 
grid used in simulations (Wang & White 2007), but even so 
they could affect the density profile of more massive haloes 
formed from their mergers and accretion. In addition, the 
expected size of the core (if any) for the relevant WDM par- 
ticle masses is close to the resolution of current simulations, 
which might explain the negative results above. 

All these uncertainties would disappear if the halo den- 
sity profile could be inferred analytically down to arbitrar- 
ily small radii directly from the power-spectrum of density 
perturbations. Salvador-Sole at al. (2012, hereafter SVMS) 
have recently built a model for the inner structure of dis- 
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sipationless collisionless dark matter haloes in bottom-up 
hierarchical cosmologief] that fills this gap. 

In the present Letter, we apply the SVMS model to 
the AWDM cosmology that results from 2 keV non-thermal 
sterile neutrinos, i.e. the minimum mass compatible with the 
Lyman-a forest (Boyarsky et al. 2009). The corresponding 
linear power-spectrum is given by (Viel et al. 2005) 



-PwDM 



(1) 



where Pcdm(&), is the linear power-spectrum for the wmap7 
concordance ACDM cosmology (Komatsu et al. 2011) and 
Tf s (fc) is the 'transfer' function, well-fitted by 



T fc (fc) = [1 + (afc) 2 ^n 5/ '\ 



with fj, = 1.12 and 

0.22 

a = 



\0.7 V 1 keV ) 



0.25 



(2) 



Mpc, (3) 



being h the current value of the Hubble parameter in units 
of 100 km s" 1 Mpc" 1 , mwDM the WDM particle mass and 
Qwdm the WDM density parameter. For the cosmology here 
considered this yields a — 0.032 Mpc. By comparing the halo 
density profile predicted in this WDM cosmology with that 
found in the underlying CDM one, we analyse separately 
the effects of the spectrum cutoff and the initial velocity 
dispersion. In Section [21 we describe the method used, in 
Section [3] to infer the density profile for halo seeds and, 
in Section [4j the density profile for haloes themselves. We 
conclude in Section [5] 



2 THE MODEL 

In any bottom-up hierarchical cosmology, haloes form 
through either major mergers or smooth accretion (includ- 
ing minor mergers). Both processes involve the virialisation 
of the halo each time the mass increases. As virialisation is 
a relaxation process yielding the memory loss of the past 
history, the density profile for haloes having suffered ma- 
jor mergers is indistinguishable from that for haloes hav- 
ing grown by pure accretion (PA; see SVMS for a complete 
rigorous explanation). Consequently, we have the right to 
concentrate in this latter kind of haloes. 

As shown in SVMS, accreting dissipationless collision- 
less dark matter haloes evolve from the inside out, keeping 
their instantaneous inner structure unaltered. In these con- 
ditions, the radius of the sphere with mass M is exactly 
equal to its virial radius, 

3GM 2 



' 10 \E P (M)\ ' (4) 

where E P (M) is the total energy of the sphere encompassing 
the same mass M in the (spherically averaged) seed, namely 
a peak in the primordial random Gaussian density field fil- 
tered at the scale M. Thus, provided the energy distribution 
in peaks is known, equation (j4| is an implicit equation for 
the halo mass profile M(r). 



1 In WDM cosmologies there is a minimum halo mass, but haloes 
still form hierarchically through the merger of less massive objects 
previously formed and their accretion of diffuse matter. 



We must remark that equation Q is only valid provided 
the isodensity contours in the seed reach turnaround with- 
out shell-crossing at increasingly larger radii (see SVMS). 
This is certainly the case when the initial peculiar velocities 
are negligible, as those induced by random Gaussian density 
fluctuations, and the seed expands in linear regime. How- 
ever, if there are peculiar velocities of non-gravitational ori- 
gin such that they dominate the dynamics at small enough 
scales, then the system may not expand in linear regime and 
there may be shell-crossing before turnaround, so equation 
may no more be valid. We will come back to this possi- 
bility in Section^ 

In the parametric form, E P (M) is given by the total 
energy in the sphere with radius R p centred at the peak, 



(•») 



E p {R p )=4n dr p r 2 <p p )(r p ) 



f [H{t^r p -v p {r p )f a 2 DM ( ti ) GM(r p ) \ 
\ 2 2 r p \ 



together with the mass of the sphere, 



M 



4tt / dr p r p (p p )(r p ) 
'o 



(6) 



where (p P )(r p ) is the spherically averaged (unconvolved) 
protohalo density profile, M(r p ) the corresponding mass 
profile, H(ti) the Hubble parameter at the cosmic time 
ti when the seed is considered, <7DM(ti) the adiabatically 
evolved particle velocity dispersion of non-gravitational ori- 
gir0 and v(r p ) is, to leading order in the deviations from 
spherical symmetry, the peculiar velocity at the radius r p in 
the seed owing to the inner mass excesq 



v p (r p ) 



2G [M(r p ) - 47rrgp(ti)/3] 
&ff(ti)rg 



(7) 



being p(ti) the mean cosmic density at U. 

The steps to be followed are thus the following ones: 1) 
determination of the spherically averaged protohalo density 
profile (p P )(r p ) from the linear power-spectrum of the cos- 
mology considered (see Sec.[3|, 2) calculation from it of the 
energy distribution E P (M) (eqs. [5]-[6]), and 3) derivation 
of the typical halo mass profile M(r) by inversion of equa- 
tion Q , and of the typical spherically averaged halo density 
profile, through the trivial relation 



(p)(r) 



1 AM 

4nr 2 dr 



(8) 



3 PROTOHALO DENSITY PROFILE 

In PA, every halo ancestor along the continuous series end- 
ing at the halo with M at f also arises from one peak in the 
random density field at ti, filtered at the mass scale of the 
ancestor. Thus, the value at the origin (r p = 0) of the con- 
volution of the spherically averaged density contrast profile 
for the protohalo, (S p )(r p ), by a Gaussian window of every 



2 The velocity dispersion due to random density fluctuations is 
several orders of magnitude smaller and can be safely neglected 
(see SVMS). 

3 In equation Q we have taken into account that the cosmic 
virial factor /(H) Ri il 0,1 is at ti very approximately equal to one. 
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Figure 1. Central density contrast of peaks in the filtered density 
field at z = 100 giving rise by PA to current haloes with masses 
equal to 10 9 M (brown lines), 10 11 M (red lines) and 10 13 M 
(orange lines) as a function of the Gaussian filtering radius Rf. 
The different lines correspond to the ACDM concordance cosmol- 
ogy (dashed lines) and the AWDM cosmology with 2 keV sterile 
neutrinos (solid lines). Radii are in physical units. 



radius R{ must be equal to the density contrast of a peak 
equally convolved, 



47T 



dr D 



2\R { ) 



(9) 



(2 7 r) 3 / 2 J R t 3 

Therefore, provided the peak trajectory, <5 p k(Rf), associated 
with the accreting halo were known, equation ([9]) could 
be seen as a Fredholm integral equation of first kind for 
(<5 p )(r p ). Such an equation can be solved (see SVMS for de- 
tails), so this would lead to the density profile (p P (r p )) for 
the seed of the halo evolving by PA. 

Furthermore, in any random Gaussian density field, the 
typical peak trajectory S p k(Rf) leading to a purely accret- 
ing halo with typical density profile is the solution of the 
differential equation (see SVMS and references therein) 

dS pk 
dR { 



= ~x e {Rt,& pk )a 2 {Ri)Ri , 



(10) 



where 02 (-Rf) is the second order spectral moment and the 
inverse of x e (R{, 5 p k) is the average inverse curvature x (mi- 
nus the Laplacian over 02) of peaks with density contrast 
<5 p k at the scale -Rf (see SVMS for the explicit form of these 
two functions). Equation (| 10f) can be solved for the bound- 
ary condition <5 p k[-Rf (M)] = 8(t) leading to the halo with 
M at t according to the one-to-one correspondence between 
peaks and haloes given by (Manrique & Salvador-Sole 1995) 



Rt(M) = 



and 



5(t) = 5 c (t) 



3M 



1 1/3 



■inp(ti) 



G(ti) 
G(t) 



(11) 



(12) 



where q is the radius, in units of Rt, of the collapsing cloud 
with volume equal to M/p(ti) associated with the peak, G(t) 
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Figure 2. Spherically averaged (unconvolved) density profile for 
the same halo seeds and cosmologies as in Figure [T] (same lines). 



is the cosmic growth factor and S c (t) is the critical linearly 
extrapolated density contrast for spherical collapse at t. In 
the underlying ACDM cosmology here considered, such a 
correspondence is given by q = 2.75 and S c (z) — 1.82 + 
(6.03 - 0.472Z + 0.0545z 2 )/(l + 0.000552z 3 ) where z is the 



redshift corresponding to t[ 

The peak trajectories, <5 p k(-Rf), solution of equation (|10|l 
in the WDM (CDM) cosmology for current haloes with three 
relevant masses are shown in Figure [T] As can be seen, the 
5 P k(-Rf) trajectories in the WDM cosmology level off, con- 
trarily to those in the CDM cosmology, at some value 5t s 
that depends on the halo mass. The time t{ s corresponding 
to Sf s (eq. [12]) marks the time when the first ancestor of 
the halo with a given current mass M forms and initiates 
the continuous series of ancestors leading by PA to that fi- 
nal halo. Before that time there is no ancestor of halos with 
that current mass in the WDM cosmology. Instead, there 
are halo ancestors down to any arbitrarily small time in the 
CDM cosmology, with no spectrum cut-off. 

In Figure [5] we show the spherically averaged density 
profile, (p P }(r p ), of WDM (and CDM) halo seeds resulting 
from equation ([9| for the peak trajectories depicted in Fig- 
ure[T] As can be seen, contrarily to their CDM counterparts, 
the WDM halo seeds show apparent flat cores with a uni- 
versal mass very approximately equal to the mass Mf s = 
[4tt/3] p(t ) (A| ff /2) 3 = 6.5 x 10 5 M associated with the ef- 
fective free-streaming scale length A|f = 2n/kf^ = a. Note 
that this is different from the mass associated with scale 
length equal to 2-7T over the wavenumber where Tf s {k) de- 
creases to 0.5, also often used to estimate the free-streaming 
mass. This latter scale length is substantially greater than a, 
so it would correspond to the scale length where the effects of 
WDM begin to be noticeable rather than the minimum size 
of density perturbations. For this reason, following Schnei- 



4 This yields the right halo mass function at every redshift (see 
Manrique & Salvador-Sole 1995). 
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der et al. (2011), the mass ~ 1.8 x 10 9 Mq associated with 
it will be called the half-mode mass, Mh m - 

If the velocity dispersion in the seed were negligible, 
its flat core with mass ~Mf s would expand and collapse 
at once, with no shell- crossing before turnaround^- Conse- 
quently, these flat protohalo cores with mass ~Mf a set a 
lower bound for the mass of WDM haloes. In contrast, there 
are no flat cores in the CDM cosmology, so there is no lower 
bound for the mass of CDM haloes either. One can find an- 
cestors of any accreting halo with M at t down to arbitrarily 
small cosmic times and with arbitrarily low masses. 



4 WDM HALO DENSITY PROFILE 

The lower-bound mass, Mf s , for haloes in the WDM cosmol- 
ogy just mentioned arises from the spectrum cutoff. But this 
bound mass has been obtained by neglecting the WDM par- 
ticle velocity dispersion, <tdm- Actually, the non-negligible 
velocity dispersion of WDM particles causes the total energy 
Ep(Rp) in the seed (eq. [5]) to become positive for radii R p 
below some value Re- This means that shells inside that ra- 
dius, with velocity dispersion larger than the Hubble veloc- 
ity, expand more rapidly than outer ones, which causes the 
system to run out of linear regime, leading to the formation 
of caustics. These caustics will fragment (owing to the per- 
turbation of the surrounding matter) and give rise to small 
nodes with masses substantially less than Me = M(Re), 
the mass of that part of the seed with null total energy. 
This means that haloes with masses below Me > Mc s do 
not actually evolve in the bottom-up fashion and that the 
typical WDM halo density profile will not satisfy equation 
@ at radii enclosing Me. 

According to Boyanovsky et al. (2008) , the value of <jdm 
today for non-thermal sterile neutrinos is related to the free- 
streaming wavenumber through 



, off 



3# 2 (t )^ 



2ct E>m(*o) 



1/2 



(13) 



In the present case, this leads to o"dm(£o) = 0.23 km s _1 . 
This value is markedly greater than the one usually adopted 
in this kind of studieqj, which might overestimate the ef- 
fects of the velocity dispersion on the predicted structure 
of WDM haloes. For this reason, we will be more conserva- 
tive. We will adopt o dm (^o) = 0.075 km s~ and study the 
effects of changing the value of ctdm(^o) by comparing the 
density profiles so obtained (WDMNT model) with those 
arising from a null value of <tdm(£o) (WDM0 model), so to 
see the effects of the cutoff in the spectrum alone, as well 
as a value of (TDM(to) equal to that of thermal neutrino- 
like WDM particles (WDMT model). According to Steffen 



5 In fact, these shells would not cross each other even after 
turnaround, meaning that such a flat perturbation would oscillate 
for ever and would not virialise. However, the subsequent shells 
coming from outside the flat core do cross them and the system 
finally virialises. 

6 Coh'n et al. 2008 take o"dm(*o) equal to 3 times the thermal 
velocity dispersion. 
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Figure 3. Typical spherically averaged density profiles predicted 
for current haloes with the same masses and in the same cosmolo- 
gies (same colours) as in the previous Figures: AWDM cosmology 
with null velocity dispersion (WDMO; solid lines), with thermal 
velocity dispersion (WDMT; short-dashed lines) and with non- 
thermal velocity dispersion (WDMNT; dotted lines); ACDM cos- 
mology (CDM long-dashed lines). To avoid crowding the profiles 
corresponding to 10 9 Mq and 10 13 Mq have been shifted 2.5 dex 
downwards and upwards, respectively. 



(2006), this latter velocity dispersion is given by, 



o-DM(io) = 0.041 3 



h \ (m 



»WDM 



keV 



"4 fi v 



0.25 



(kr 



(14) 



leading for 2 keV particles to (tdm(^o) = 0.015 km s _1 . 
The reference ACDM halo profiles are derived assuming, of 
course, a negligible ctdm(^o). 

In Figure [3] we plot the typical spherically averaged 
halo density profiles for current haloes with the same masses 
as used in the previous Figures, each of them for the three 
values of (Xdm(£o) just mentioned. All the profiles devi- 
ate from the corresponding CDM halo profiles in the same 
monotonous way leading to a flat core with mass Mf s . How- 
ever, the only profile that can be strictly traced down to 
the halo centre is the WDMO profile. The remaining pro- 
files (for non- vanishing odm^o)) show a sharp up-turn to 
infinity at a small enough radius. This reflects the fact that, 
for M approaching the value Me where E P (M) vanishes 
(eq. [5]), equation Q is no longer valid. For thermal veloc- 
ity dispersion (WDMT profile), the minimum radius reached 
essentially coincides with the edge of the flat core, while, for 
larger values of crDM(to), it is significantly larger. The mass 
Me encompassed by the minimum radius reached in the case 
of non-thermal velocity dispersion (WDMNT profile) , of the 
order of the half-mode mass, Mhm, for the value of UDM(to) 
used here, shows a slight trend to diminish with increasing 
halo mass. Specifically, it is equal to 10 9 Mq, 1.9 x 10 8 M 



and 8.1 x 10 7 M Q for haloes with 10 9 M G 
Mq, respectively. 



10 11 M and 10 13 
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5 CONCLUSIONS 

The typical spherically averaged density profile for haloes 
with masses greater than the free-streaming mass, Mf s , in 
the WDM cosmology can be derived analytically by means of 
the SVMS model. As a consequence of the spectrum cutoff, 
the density profiles so obtained show a clear flattening rela- 
tive to the profiles found in CDM cosmologies, independent 
of the particle velocity dispersion, which gradually evolves 
into a flat core with mass equal to Mf B . A particle mass as 
low as 2 keV is enough for the flattening to begin to be vis- 
ible, at a radius of ~ 0.005 virial radii as reached in current 
./V-body simulations, for haloes with masses M as large as 
~ 10 13 M and the flat core is already reached, at those 
radii, for haloes with M < 10 11 M . 

The only apparent effect of the particle velocity dis- 
persion is that it prevents from reaching the flat core. In 
fact, such a core might not exist, that is the density pro- 
file at the halo centre might not be flat. Indeed, the veloc- 
ity dispersion causes protohaloes to have null total energy 
within some small radius encompassing the mass Mf s . This 
produces caustics as protohaloes expand, leading to their 
fragmentation into small nodes. This would explain the pres- 
ence of haloes with masses below Mf s in TV-body simulations 
of WDM cosmologies. Therefore, small mass haloes do not 
grow in the bottom-up hierarchical fashion and their density 
profile cannot be recovered by means of the SVMS model. 
For a 2 keV particle mass and a non-thermal velocity dis- 
persion of ~ 0.075 km s _1 , the "unresolved core" affects the 
whole density profile for haloes with masses below ~ 10 9 
Mq. The larger the velocity dispersion, the larger the mini- 
mum radius that can be reached in the density profile of any 
given halo. All these effects are obviously relaxed for increas- 
ing particle masses, which cannot be ruled out according to 
the present results. 
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